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1 Overview 


This report summarizes the activities on NASA Grant NAG- 1-2249 during the past year. 

The work summarized below consists of two parts. The first part concentrates on the lumped 

element model of synthetic jets, while the second part focuses on the structural dynamic 

modeling of piezoceramic flap actuators. 

1.1 Objectives 

• Develop coupled electro/fluid/structural lumped-element model (LEM) of a prototypical 
flow-control actuator 

• Validate the coupled electro/fluid/stmctural dynamics lumped-element models 

• Develop simple, yet effective, design tools for actuators 

• Develop structural dynamic models that accurately characterize the dynamic response of 
piezoelectric flap actuators using the Finite Element Method (FEM) as well as analytical 
methods 

• Perform a parametric study of a piezo-composite flap actuator. 

• Develop an optimization scheme for maximizing the actuator performance. 

1.2 Progress 

• Developed electro/fluid/structural models of piezoelectric-driven synthetic jet actuators using 
lumped-elements modeling 

• Designed and fabricated a modular setup of a synthetic jet actuator driven by a piezoelectric- 
diaphragm to conduct careful experiments to isolate model components 

• Designed optimized piezoelectric-diaphragms 
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• Developed structural dynamic model based on FEM for the design of piezoelectric flap 
actuators 

• Designed, fabricated, and characterized the flap actuators 

• Developed an analytical scaling method and an optimization scheme that enables to select the 
actuator parameters to maximize its performance for a given bandwidth. 

2 Synthetic Jet Actuators 

2.1 Theoretical Model of a Piezoelectric-Driven Synthetic Jet Actuator 

A coupled electro-mechanical-acoustic model of a synthetic jet with piezoelectric- 

diaphragms was developed. The main assumption employed in lumped-element modeling 
(LEM) is that the characteristic length scale(s) of the governing physical phenomena are much 
larger than the largest geometric dimension. The governing equations for the distributed system 
are “lumped” into a set of coupled differential equations. Device components are represented by 
an equivalent circuit element and are coimected by following specific rules based on the effort or 
flow variable type they share. The benefit of the LEM lies in the explicit relationship between 
each device component and the frequency response of the system. Such an approach provides a 
simple method to estimate the dynamic response of a synthetic jet and ultimately to design it. 

The actuating diaphragm consists of a clamped axisymmetric piezoelectric unimorph actuator 
in which a piezoceramic is bonded to the central portion of a metal diaphragm. The composite 
diaphragm is driven into motion via an applied AC voltage. This represents a conversion from 
the electric to the mechanical domain. Since the primary purpose of the diaphragm is to produce 
large volume displacements, we focus on the first axisymmetric vibration mode of the composite 
diaphragm. To accomplish this, the diaphragm is lumped into an equivalent mass representing 
the stored kinetic energy and a compliance representing the stored potential energy. The 
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effective mechanical mass and compliance can be easily converted to their acoustic coimterparts, 
and so it is convenient to treat the entire device as a couple electroacoustic system. Since the 
cavity contains a compressible gas, it too stores potential energy and it is modeled as an acoustic 
compliance. Finally, the orifice is a viscous region associated with motion of a fluid mass. 
Therefore, there will be an effective mass and a resistance associated with the orifice neck. The 
neck acoustic resistance takes into accoimt the viscous loss in the neck and the non-linear 
discharge coefficient of the jet exit that is approximated by modeling the orifice as a generalized 
Bernoulli obstraction. For a first assumption, fully developed laminar pipe flow with oscillating 
pressure gradient is investigated to model the flow in the orifice. In addition, for low operational 
frequencies, a radiation mass is modeled, to first order, as a piston in an infinite baffle if the 
circular orifice is flush-mounted in a plate that is much larger in extent than the orifice radius. 

The equivalent circuit representation of a synthetic jet is shown below. Mom and CaM are the 
acoustic mass and compliance of the piezo-diaphragm, respectively. Com is the acoustic 
compliance of the cavity, while Rao and Mao are the acoustic resistance and mass of the orifice, 
respectively. and MaRad are the acoustic resistance and mass associated with the orifice 

radiation, respectively. The conversion fi:om the electrical to the acoustic domains is accoimted 
for via a transformer with a turns ratio <^a- This converts the driving current to a volume velocity 
of the diaphragm. The goal of the design is to maximize Qm/Vac (i e-, the volume flow rate 
through the orifice per applied voltage). 
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Figure 1: Acoustical equivalent circuit representation of a piezoelectric-driven synthetic jet 
The transfer function of the whole equivalent circuit is given by 





a^s* + + 0,5 + 1 


, where 


Oj = {R^ +^a/) 

a, = C^yM^C^ M.h ) 

Although this expression appears quite complex, it reveals some important features of the 
synthetic jet. First, the frequency response function is the ratio of volume velocity to the applied 
voltage. Since volume velocity is the derivative of the volume displacement, we see that for low 
frequencies (s —> 0), the volmne velocity is zero and the volume displacement is (j>a CaM=-da, as 
expected. Second, the denominator is a 4*-order polynomial in s, indicating two resonance 
frequencies. It can be shown that these two resonance frequencies are controlled by the 
diaphragm natural frequency fd and the Helmholtz frequency of the cavity fc, given by 


such that 




2;r 


/c= — 


2;rV(A/^+M^)C, 


( 2 ) 
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( 3 ) 


2.2 Design of Modular Synthetic Jet Actuator 

A modular apparatus of a piezoelectric-driven synthetic jet with cylindrical orifice was 

designed and constructed in order to perform careful experiments of each device component. 
The modular design allows for the interchanging of parts to test a variety of geometries, 
including the orifice, cavity, and different piezoelectric-diaphragms. Two mounts accommodate 
different composite piezoelectric-diaphragms (diameter of 27 mm or 41 mm, along with various 
thickness for both the piezoceramic and the shim). In addition a dual configuration can be 
obtained that is, two piezo-diaphragms are mounted face-to-face. The size of the cavity is 
modular in the same way with its height varying from 2 to 6 mm. Finally, various orifice plates 
provide an orifice with a diameter ranging from 0.28 mm to 0.78 mm and a depth fi'om 0.5 to 2 
mm. To mormt the microphone flush against the bottom of the body plate, a tapered hole was 
machined from the end of the mount up to the cavity, which permitted insertion of the 
microphone without allowing air to escape. The basic device is shown below in Figure 2, along 
with one picture in Figure 3. 

OrificepUte 


Top pUte Body pUte 

Access hole 
for 

. xnicrophoBe 
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Figure 2: Schematic of the modular synthetic jet apparatus. 



Figure 3: Picture of the synthetic jet tested. 


Figure 4; NASA synthetic jet device, courtesy 
of Dr. Michael Kegerise. 


In addition to this modular design. Dr. Michael Kegerise of NASA Langley provided us with 
the design drawing of a synthetic jet device that is currently under investigation, as shown above 
in Figure 4. This one differs essentially from the previous device by its rectangular slot, 
compared to the cylindrical orifice. 

2.3 Experimental Procedures 

In order to verify each component of the LEM, the experimental apparatus built permits a 
systematic parametric study of the synthetic jet design. The procedures are as follows: 

• Piezoelectric-diaphragm excitation in a vacuum chamber to isolate the 
electromechanical coupling 

• Oscillatory pressure in a sealed cavity to eliminate orifice effects 

• Behavior of synthesized jets versus various orifice geometries to determine the orifice 
impedance in ambient flow conditions 
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2.4 Experimental Results 

To gain confidence in the use of LEM as a potential tool to optimize the design of synthetic 
jets, we applied the above LEM approach to the NASA synthetic jet design since preliminary 
performance data was already available for that device. The predicted fi-equency response 
fimction between the applied voltage and the output volume flow rate is shown in Figure 5. The 
predicted resonant frequencies are 435 and 1557 Hz, while the estimated natural fi-equencies of 
the cavity and diaphragm are 443 and 1527 Hz. 

Figure 6 shows the results of hot-wire measurements conducted by NASA to characterize the 
maximum velocity. Two curves are shown, one each for the maximum RMS and maximum 
amplitude of the velocity measured near the slot. The agreement between the LEM and the 
experiment is quite promising. It can be noted that, due to the lumped nature of the LEM, only 
the average velocity can be estimated. Clearly, the magnitude of the response at the resonance 
fi-equencies is not predicted correctly. This is undoubtedly due to the oversimplification of the 
resistance terms in Jiao- Nonetheless, the LEM predicts, using a driving AC voltage of 1 12 V, an 
average velocity in the range of 800-1000 Hz of approximately 8 m/s. This is in good agreement 
with experiment. 
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Also, experiments of the piezoelectric-diaphragm itself were conducted to validate the 
composite linear plate theory used to model it. A maximum input voltage of 5 V drove the 
piezoelectric-diaphragm into motion, at a frequency on the order of 100-200 Hz, and the 
deflection of the di^hragm were measured by a scanning laser vibrometer. Figure 7 below 
shows the comparison between the experiment and the dieory. 



Figure 7: Volume displacement of the piezoelectric-diaphragm by applied voltage. 

Similarly, to test the theory on the cavity acoustic compliance, the orifice was replaced with a 
solid cap to provide a closed cavity and all leaks were carefully minimized. The piezoceramic 
was then driven with a nominal 1 V amplitude sinusoid, and the displacement of the vibrating 
diaphragm was measured with a laser displacement sensor (Micro-Epsilon Model ILD2000-10). 
A 1/8 in. Briiel & Kjaer (B&K) type 4138 condenser microphone with B&K type 2669 
preamplifier measured the fluctuating pressure in the closed cavity. The amplitude of the 
sinusoid was adjusted to avoid harmonics. The frequency response function between die 
pressure and displacement signal at the diaphragm center was used to measure P/w^ at 
several frequencies and calculate C^y, as shown in Figure 8. Using the average measured value 


of 9.89 AfPa/m and the measured mode shape, the cavity volume was determined to be 
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2.53x10-6 ± 12%. This nominal value is within ~1% of the cavity volume calculated from 

the geometry. 
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Figure 8: Measured acoustic compliance C^y vs. frequency in closed cavity of synthetic jet 
Finally, the lumped element model is used to predict the frequency response of two 
prototypical synthetic jet actuators. The dimensions and material properties of the piezoceramic 
diaphragms and actuator devices, along with the pertinent lumped element parameters, can be 
found in details in Gallas et al. (2003). Figure 9 show the model prediction of the centerline 
velocity compared to phase-locked LDV measurements versus frequency. 



Figure 9: Comparison between the lumped element model and experimental frequency response measured 
using phase-locked LDV for two prototypical synthetic jets (Gallas et aL 2003). 

2.5 Future Work 

Future work will focus on the validation of each device component. Furthermore, we need to 
look in the acoustic liner literature for model improvement of the orifice behavior, including 
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nonlinear and grazing flow effects. In the same way, experiments and LEM results will be 
comparing with ongoing CFD simulations of the orifice flow conducted by Dr. Raj at Mittal of 
GWU. Once the LEM approach is validated, it will be used as basis for a synthetic jet design 
tool for end users utilization. 

2.6 Publications/Presentations 

Gallas, Q., Holman, R., Nishida, T., Carroll, B., Sheplak, M., and Cattafesta, L., “Lumped 

Element Modeling of Piezoelectric-Driven Synthetic Jet KcXxxztots” AIAA Journal, No\. 
41, No. 2, pp. 240-247, February 2003. 

Gallas, Q., Mathew, J., Kaysap, A., Nishida, T., Sheplak, M., and Cattafesta, L., “Lumped 

Element Modeling of Piezoelectric-Driven Synthetic Jet Actuators, AIAA Paper 2002- 
0125, 40th AIAA Aerospace Sciences Meeting &. Exhibit, January 2002. 

Gallas, Q., Sheplak, M., and Cattafesta, L., “Lumped Element Modeling of Synthetic Jets, Ultra- 
Efficient Engine Technology Forum, NASA Gleim Research Center, September 2001. 

3 Piezoelectric Flap Actuators 

Piezoelectric flap actuators are desirable for active flow control applications. An 
understanding of their characteristics is essential prior to using them in a fluid flow. This portion 
of the report summarizes the theoretical modeling and experimental validation of the dynamic 
response of piezoelectric unimorph actuators. Two different finite element models (FEM) of a 
piezoelectric unimorph actuator are developed. One is a beam model that assumes a perfect 
bond exists between the piezoelectric patch and the shim, and the second extends the perfect 
bond model by incorporating a linear elastic shear element for the bond layer. These models are 
then used to predict the magnitude of the dc response (i.e., tip displacement per imit applied 
voltage) and the natural fi-equency of a piezoelectric flap. Also an approximate analytical static 


13 



model was developed to determine the functional dependence of the system response on the 
actuator parameters. The analytical model was also extended to study the response of bimorph 
actuators. The results are compared with experimental data obtained jfrom a parametric study, in 
which ten otherwise identical piezoelectric flaps with varying piezoelectric patch sizes are 
fabricated and characterized using a laser displacement sensor. The results indicate that the 
theoretical models provide good estimates for the dc response and for the natural frequency when 
a linear elastic bond layer is used and a clamped boundary condition is realized. In particular all 
models produce estimates that are accurate to within 15% and 10% for die dc response and 
natural frequency, respectively. Finally a design procedure for developing actuators with 
optimum performance is briefly summarized. 

3.1 Theoretical Model of a Piezoelectric-Driven Flap Actuator 

A parametric study was carried out in which the geometry of a piezoelectric-shim was varied 

from approximately 10% to 100% of the total beam length in ten equal steps (while all other 
variables were fixed). The flaps were constructed of an aluminum alloy shim, with the 
piezoelectric material firmly affixed to its surface using an epoxy adhesive to form a piezo- 
composite beam. Figure 3 shows the schematic of a piezoelectric flap unimorph actuator. The 
total length of the shim is Lj. (~70 mm), the width is b (~26 mm), and the thickness is (508 

pm). A piezoceramic of length (variable) and thickness (254 pm) is bonded to the shim 

using an epoxy adhesive of thickness (13-25 pm). A small gap of length (< 1 mm) is 

maintained between the clamp and the piezoelectric to eliminate stress concentrations and to 
prevent an electrical short circuit. Upon application of an external voltage to the piezoelectric 
electrode, the flap deflects in a direction determined by the orientation of the piezoelectric 
polarization vector with respect to the electric field. 
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In order to characterize piezoelectric flap actuators, a structural dynamics model is required 
that can accurately predict the frequency response of the flaps. Since flow-control applications 
of interest typically demand a large volume displacement, attention is restricted to the first 
bending mode only, which corresponds to the maximum volume displacement. Below, a 
“perfect bond” model and a “shear lag” model, both of which utilize finite element 
methodologies are used to calculate the frequency response of the actuators to an applied 
voltage. A simpler analytical static analysis is then described that is suitable for design. 

The quantities of interest, namely the dc response and natural frequency of the actuator were 
calculated using the FEM and scaling analysis, and compared with the experimental results. 


z 



Figure 3: Flap actuator model. 

3.2 Perfect Bond Finite Element Model 

A simple schematic of the perfect bond model is shown in Figure 4. This assumes an 
infinitely rigid bond exists between the piezoelectric and the shim. All loads applied to the bond 
by the piezoelectric are transmitted fully to the shim. 
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Crawley and de Luis [1987] showed that a perfect bond between the piezoelectric and the 
shim leads to a concentrated couple or moment at the edges of the piezoelectric patch 

X = and . As described in Cattafesta et al. [1997], a relationship can be obtained 

between the applied voltage and , assuming that the actuator is in pure bending, and is given 
by 


^o~ 2 (^‘'2 H) 

where is the elastic modulus of the piezoelectric, J 3 , is the appropriate piezoelectric constant 

for a bender configuration, is the magnitude of the applied electric field, b is the width of 

the flap, tp is the piezoelectric thickness, and is the location of the neutral axis in the 

composite section with respect to the bottom surface of the piezoelectric and is given by 


c, = — 
^ 2 




1 .^ 
tj_ 

Et, 


(5) 


where E^ is the elastic modulus of the shim and is the shim thickness. 
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Figure 4: Perfect bond model. 

It is important to note that a perfect bond cannot be achieved in practice because of the finite 
stiffiiess of the bond. Hence, a shear stress exists at the bond interface. This situation is 
considered in the next section. 

A standard Euler-Bemoulli beam finite element model was developed to calculate the grid- 
resolved frequency response fiinction for the piezoelectric actuators. The details of the finite 
element modeling scheme are given in Cattafesta et al. [2001]. The material properties of the 
flap are summarized in Table 1. 
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Table 1: Flap Actuator Properties. 


Property 

AL-2024 

(Aluminum 

shim) 

PZT-5H 

(piezoceramic) 

Eccobond 45 
(Adhesive) 

Elastic Modulus 
(Pa) 

7.33T0*® 

6.2T0‘° 

3.9810® 

Density 

(kg/m^) 

3290 

7500 

1060 

Piezoelectric 
constant (mAO 

N/A 

-274T0'*^ 

N/A 

Shear Modulus 
(Pa) 

N/A 

N/A 

1.4810® 


3.3 Shear Lag Finite Element Model 

The perfect bond model described above results in an artificially stiff model for the flap 
actuator, causing predicted values for the natural fi-equency that are higher than the 
corresponding measured values (Cattafesta et al. [2001]). A more realistic model for the bond is 
thus sought to obtain better quantitative predictions for the flap characteristics. 

A perfect bond model neglects the presence of the bonding adhesive layer between die 
aluminum shim and the piezoelectric. The bond layer is viscoelastic to some extent and 
possesses a finite stiffiiess. Some amount of strain will occur in the bond layer, resulting in a 
shear lag between the induced strain by the piezoelectric and the strain actually transmitted to the 
shim via the bond layer. 


18 





















Crawley and de Luis [1987] showed that the maximum shear stress occurs near the two edges 
of the piezoceramic patch, whereas the shear stress distribution is almost zero for the remaining 
regions of the bond. The stress distribution also has opposite signs at the two edges of the bond. 
The shear-lag model described below represents an extension of their work to investigate the 
effect of the bond layer on the dynamic response of the actuators. 

A schematic of the shear lag model is shown in Figure 5. Due to the deformation of the 
piezoelectric induced by the application of an external electric field, a horizontal 
electromechanical force of magnitude / = acts on the front and back faces of the 


piezoelectric section, where E^ is the elastic modulus and A = bL^ is the cross-sectional area of 


the piezoelectric. A shear stress t(x) acts at the interface between the bond and the 
piezoelectric (and the bond and shim), as shown in Figure 5. A fi'ee-body diagram of the 
piezoelectric reveals the relationship between the magnitudes of / and r(x). 


r /2 

’’ r{x)bdx . An inspection of Figure 4 for the perfect bond case reveals the effect of this 


force is to produce an induced couple in Eq. (4). The relation between the horizontal 
electromechanical force and the induced couple can be obtained by combining Eq. (4) with the 
definition for / and is given by the formula 


/ = 




( 6 ) 


Thus, the flap deflects or bends upon application of a voltage, and the direction of motion is 
determined by the orientation of the piezoelectric polarization vector with respect to the electric 
field. The net effect is that a finite shear stress or shear lag, r(x) = G{dujdz) , is induced in the 


bond layer due to the deformation of piezoelectric, resulting in a decreased induced strain being 
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transmitted to the shim. Here, u is the displacement in the x direction and G is the rigidity 
modulus. 

A finite element model, similar to the one used for the perfect bond model, was developed. 
The solution procedure is same as that for the perfect bond model, the details of which are given 
in Mathew et al. [2001]. 



Figures: Shear lag modeL 

3.4 Analytical Scaling Analysis 

z 
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In order to obtain the functional dependence of the critical quantities like dc response and 
natural frequency on the flap design parameters, an analytical scaling analysis was undertaken. 
The sections below briefly summarize the modeling scheme used. 


3.4.1 Calculation of dc response 

The Euler-Bemoulli beam equation for the flap was solved for an applied external voltage. 
Figure 6 shows the arrangement of a unimorph piezoelectric patch that is perfectly bonded to a 
shim material. The application of a voltage produces a couple that acts at the edges of the 
piezoelectric patch, as described in the perfect bond model. The composite beam is split into 
three sections. Section 1, of length is the small gap region between the clamp and the 

piezoelectric. Section 2, of length , is the composite section. Section 3, of length , is the 
second shim section. 

In Section 1, the governing static equilibrium equation for the vertical deflection of the beam 
is 


^ dx^ ' 


In Section 2, the governing static equilibrium equation is given by 


(7) 


dx^ 


+EJ, 


u vv. 


^ = M 


( 8 ) 


where is due to the applied voltage and is given by Eq. (4). In Section 3, the governing 
static equilibrium equation is given by 




(9) 


To find the beam deflection as a result of an applied electric field, we use the following four 
boundary conditions for a cantilever beam 
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Wy{X = 0) = 0, 

w[(x = 0) = 0, 

^^^,(x = L,) = 0, 

and 

<(x = Z^) = 0. (10) 

These boundary conditions states that the deflection and slope are zero at the clamped edge and 

the force and moment are zero at the free edge. The matching condition matches the deflection 

and slope at the interface between sections 


^[(.x = L^) = w[{x = L^^), 
= W3(x = 


and 


W2(^ = -^*,p+ip) = ^(^ = ^g<v+^p)- 01) 

In the ideal case the gap region is negligible in size and neglected, -> 0, so we can write 

the piecewise continuous beam deflection equation as 




w, 




Q<x^L. 


(Ell 2 ■ 


gap 




and 








(^4 


X — 




(12) 


To obtain the tip deflection, we substitute iox x = Lp, (£7)^ expression 

for Mj from Eq. (4) into the expression for Wj(x) 
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( 13 ) 


(£,V+4£a’-S/,+6£,(/£,</ + 4£,V£,(,+V,‘)’ 

It can be shown that the dc response depends on non-dimensional parameters like 
Ep!E^,tJLj,zxA d^y It^ , where Ej^ = F/r^ . Define T = tp!t^, 1 = 1^1 Lj. Tmd E = E^ IE, 

for the purposes of scaling the dc response. Substituting these into Eq. (13) and dividing by 
d^yiyltp gives 


-3£(l + F)<(L-2)i^ 

^ ^ = —L (14) 

d^yij. (l + 4Et +6Et^+ 4t ^E + E^V ) 

The above expression gives the tip deflection for a given piezoelectric material of a specified 
total length for an applied electric field. The implications of the above expression will be 
discussed at a later stage in this report. 

A similar analysis as above can be done for a symmetric bimorph actuator. A bimorph 
actuator has oppositely poled piezoelectric patches attached to both sides of the shim and is used 
in applications that require greater volume displacements. The schematic of a bimorph actuator 
is shown in Figure 7. The dimensioning convention is same as that for the unimorph actuator. 



i 
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Figure 7: Schematic of a bimorph actuator. 

The induced moment at the piezoelectric s him interface is given by the general expression 
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^0 = \y<^{yW 05 ) 

where y is the coordinate along the cross section measured from the neutral axis, cr{y) is the 

normal stress distribution, and dA is the area element for integration. Eq. (15) can be simplified 
to 


r 


^0 ~ 


t. 


( 16 ) 


Note that the neutral axis passes along the centroid of the cross-section. Following the same 


procedure as that for the unimorph, the non-dimensional tip deflection can be expressed as 


W = 




-6E{\+T)T[L-1)L^ 



[U6ET + %Ve)‘ 


( 17 ) 


3.4.2 Calculation of natural frequency 

The natural frequency of the composite beam is calculated based on the static deflection of 
the composite beam imder self load (i.e., weight). For simplicity, the gap region is neglected in 
this analysis. The schematic of the beam is given in Figure 8. In the figure, and are the 


load per unit length of the composite section and the shim section respectively. 





Figure 8: Schematic of the beam under self load. 

A free body analysis of the beam gives the force and moment acting at the clamp point as 


24 


(18) 


^ = <hL,-^Qi(Lr-L,) 


M = - 


2)^ 


, g.(V+ V) 


(19) 


Using the same nomenclature described in the previous analysis, the governing equation for 
vertical deflection in Section 2 is given as 


(et +ei\ 

\ s S2 P Pi J 




dx' 


-=M + Rx-q2 — . 

' 2 


In Section 3, the governing static equilibrium equation is given by 


( 20 ) 


El 




(x^ 


( 21 ) 


To find the beam deflection as a result of an applied electric field, we use the same boundary 
conditions, Eq. (10), and matching conditions, Eq. (11). This result in the following piecewise 
continuous beam deflection relation: 


^2 (^) = 24(^7) + ) (22) 


and 


where 



-1 

24{Ell 


^Ax"* +Bx^ + Cx^ +Dx+E^, 


(23) 


{ei\ 

D = ^(‘W. V + 12i,92 V - 12 V«>V ) + 12i,«A (i, - -tr ) - 4?, V P4) 

E = V ^ (8«2iA - 2?z (i/ +3 V ) - ) + «. V (SV -8i, A- + ) 

The tip deflection for an applied force is given by the expression 
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(25) 


t. 


^^(V H- V -2iA(2V -3i,L, +2i/)) 

As a check, note that Eq. (25) reduces to -3^2 A-"* /(24(E'/)^ ) , which is the tip deflection of a 

cantilever beam imder self load, for the limiting cases when (SI)^ = and L^=0. The 
potential energy associated with the deflection in the beam is given by the expression 


1 Zt 

U = - \E(x)I(x) 


dx^ 


dx. 


(26) 


Similarly, the kinetic energy is given by the expression 


1 ^ 2 

r=- Jm,(x)(w(x)) dx. (27) 

2 0 

where m,{x) is the mass per unit length of the actuator. Integrating Eq. (26) and Eq. (27) using 
the piecewise continuous solution, Eq. (22)- Eq. (23), we obtain the total potential and kinetic 
energy of the beam. Lumping the potential energy and the kinetic energy at the tip, we can 
obtain the effective compliance and mass of the composite beam. 

The effective compliance of the composite beam is 



(28) 


Similarly, the effective mass for the beam is obtained by the expression 



The first resonant frequency of the composite beam is estimated by 


(29) 
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The natural frequency is estimated for both the unimorph and bimorph actuators. The 
expression for the natural frequency is very long and is therefore not listed here. The natural 
frequency can be normalized by the natural frequency of a homogenous beam under self load 
and is given by the formula 


L 

IaV 

0.16221 

1 EJs" ' 


It can be shown that the non-dimensional natural frequency is a function of non-dimensional 
parameters like L^tL^, £/£^,and p^l p^. 

3.5 Experimental Setup and Procedure 

Experiments were conducted to characterize the piezoelectric actuators and obtain their 
natural frequency and dc response. The response of the actuators to an applied input voltage was 
measured. A laser displacement sensor was used to obtain the frequency response function 
between input voltage and output tip displacement. 

The set of actuators designed had a typical length of 70 mm and a width of 26 mm. The shim 
was 508 pm thick, while the piezoelectric was 254 pm thick. The piezoelectric patch was firmly 
bonded to the shim using Eccobond 45 epoxy adhesive. Eccobond 45 epoxy has a glass 
transition temperature of 48° and is not very viscoeleastic leading to minimal dissipation of 
energy as heat. The thickness of the bond layer varied slightly from 13-25 pm over the bonding 
area. Ten of these actuators were built, where the length of the piezoelectric patch was increased 
in steps of 10% between the successive actuators. The actuators were clamped at one end on a 
test apparatus to approximate a clamped-free boundary condition. Care was taken to ensure the 
rigidity of the clamp so that there is minimal motion at the base of the clamp when the actuator is 


vibrating. The actuators were mounted vertically on the test apparatus and the electric lead 
connections were realized via low resistance (<10 ohms) strain-relieved, conductive copper tape 
attached to the piezoelectric patch as well as to the grounded apparatus. 

A periodic sweep excitation signal was generated using a function generator (HP model 
E1441 A). The sweep signal spanned a frequency range from 0.1 to 200 Hz with amplitude of 
0.05 V. The voltage signal was amplified by a nominal gain of 30 using a piezoelectric amplifier 
(PCB Model 790A06). The period of the chirp was set at 4 seconds to provide good frequency 
resolution ( A/" = 0.25 Hz) for the spectral analysis. The output data was acquired using a 16-bit 
data acquisition system (HP model E1433A). Fifty block averages of the spectral data were 
formd suitable to reduce the random error component to negligible levels. The laser beam was 
focused at a point approximately 1 mm from the midpoint of the tip of the flap. The schematic 
of the experimental setup is shown in Figure 9. 



Figure 9: Schematic of the experimental setup. 
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3.6 Results and Discussion 

A laser displacement sensor was used to characterize the actuators. The frequency response 
function of the actuators was measured and the natural frequency and dc response was calculated 
from the frequency response. The results were compared with the values obtained from the two 
finite element models as well as the analytical model described previously. The dc response was 
calculated from a 2”‘*-order system fit as well as from an average of first ten frequency response 
function magnitude values. The mathematical expression for the fit is given as 


H{a>) = 


X{a>) 


1 


( 32 ) 


V -Mco^ + Cja + K 

where the frequency response function fit, H(o}) , is defined as the ratio of the tip displacement 
X(co) to the applied voltage V , for each frequency eo . The coefficients M , C, and K are the 
effective mass, damping and stiffiiess for the system respectively. The natural frequency was 
calculated from a 90° phase angle, the frequency corresponding to the peak of the power spectral 
density, theoretical fit, and an impact hammer test. The response of the actuators was calculated 
using multiple methods to obtain a good agreement between experiment and theoretical values. 
The details of these multiple methods for the calculation of dc response and natural frequency 
are given in Mathew [2002]. 

A comparison of the natural frequency for the flaps between experiment and theoretical 
values is plotted versus normalized piezoelectric length in Figure 10. The experimental results 
agree well with the theory to within +/- 10%. Theoretical values include predictions from the 
analytical model, perfect bond finite element model analysis and shear-lag finite element model 
analysis. The prediction from the shear lag model is closest to the experimental values. The 
perfect bond frequency predictions are higher than the shear lag values, while the static analysis 
data serves as an upper limit for the theoretical predictions. Note that analytical model gives 
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identical results to the perfect bond finite element model if 0 is not assumed. A peak in the 
system bandwidth is observed for = 0.5. For low values ofl^/Zj. , increasing L^jLj. has 

a greater impact on the stiffiiess than on the mass of the actuator. However, at higher values of 
LpjLj. , increasing LpjLj. has a greater impact on the mass than on the stiffiiess of the actuator. 

Thus, there is some length L^jLj. that balances these two competing effects. This gives an idea 

of the optimum piezoelectric patch length for maximizing the natural frequency of the present 
piezoelectric actuator. 

Figure 11 shows the dc response comparison between experiment and theoretical models. 
All theoretical models predict the dc response to a reasonable accuracy of +/- 15%. The dc 
response increases monotonically with increasing piezoelectric length and saturates around a 
value of LpjLj. =0.9. It is clear that increasing the piezoelectric length induces more strain, 

which leads to larger tip displacements and hence a larger dc response. It can also be observed 
that the dc response predicted by the perfect bond finite element model nearly matches that 
predicted by the analytical model. Thus the experimental data presented here give confidence in 
the theoretical predictions by the models developed, with the shear lag finite element model 
being the most accurate, while the analytical model is the simplest. 

The analytical and finite element formulation treats the actuators as a beam model. Actuator 
dimensions resemble that of a plate, which has less stiffiiess than a beam. Leissa [1993] defined 

a frequency parameter given by , where a> is the natural frequency in radians, / 

is the total length, p is the density, E is the elastic modulus and t the thickness of a 
homogeneous cantilever plate. Approximating the actuator designs by homogeneous plates, and 
calculating the frequency parameters yielded values of 3.63. However the frequency parameter 
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value for a homogeneous cantilever beam is 3.52 (Thomson [1993]). Since the frequency 
parameter values are within 3 %, approximating the actuator design by a beam model is justified. 



Figure 10: Comparison of natural frequency between experiment and theory. 


X 10"* 



Figure 11: Comparison of dc response between experiment and theory. 
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3.7 Actuator Design Considerations 

It can be observed from the experimental results that the dc response increases monotonically 
with the piezoelectric patch length and there is an optimum piezoelectric length that maximiz es 
the system bandwidth. The study revealed the dependence of the dc response and the natural 
frequency of the actuator on only one parameter, the ratio of the piezoelectric patch length to the 
total length of the actuator, L^jL^ . However, the analytical static analysis, described earlier, 

showed that the nondimensional dc response depends on non-dimensional parameters like 
LplLj., EpIE^, tJLj and the nondimensional natural frequency is a function of non- 
dimensional parameters like tpits, L^/L^, £p/^,,and 

A nominal design with L^l Lj.=0.5, E^l E^=Q.^5 , =0.0073, and 

Pp! Ps- 2-48 is selected to study the dependence of both non-dimensional dc response, W , and 

natural frequency, F , on the actuator parameters. Figure 12 (a-d) shows the variation of the 
non-dimensional dc response. For each plot only the parameter along the x axis is varied while 
all other parameters are fixed at their nominal value. 

W increases monotonically with the stif&ess of the piezoelectric patch as well as its length, 
and saturates for a value of L^l Lj= 0.9 . However, W increases initially as the thickness of the 

piezoelectric is increased until a value of 0.54, and then decreases on further increasing 

the thickness ratio. As tjL^. ratio is increased, W drops. The same behavior is shown by 
bimorph actuator in Figure 13 (a-d), only difference being that W peaks at = 0.39 . Also as 
shown in plots 11c) and lid), W drops with increasing tjL^. and increases with increasing 
!Lp , until it saturates around a value of X/I, = 0.9 . 
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The variation of non-dimensional natural frequency with the design parameters for a 
uiumorph is shown in Figure 14 (a-d). It can be seen that F increases with increasing stif&iess 
ratio and thickness ratio but drops witii increasing density ratio. F peaks at a length ratio of 
Lp! Lj. =0.53 and then drops on further increasing the piezoelectric patch length. The same 
behavior is exhibited by bimorph in Figure 15 (a-d). Here F peaks for a length ratio of 
LpIL^ =0.55. 






Figure 12 (a-d): Variation of non-dimensional dc response with the design parameters for a unimorph. 
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Figure 15 (a-d): Variation of non-dimensional natural frequency with the design parameters for a bimorph. 

It is clear from the above plots that since the actuator performance depends on a number of 
design parameters, a design optimization has to be performed in order to develop an actuator that 
gives the best performance. We are specifically interested in designing an actuator that gives the 
maximum tip deflection for a given system bandwidth, which is approximated by the natural 
frequency of the actuator. Consider the case of one of the rmimorph actuators developed that has 
a total length of 70 mm, piezoelectric patch length of 35 mm Lj.= 0.5), shim thickness of 

508 pm, and piezoelectric thickness of 254 pm. The natural frequency of this actuator calculated 
from analytical scaling analysis is 137 Hz and its non-dimensional dc response, W is 44.2. For 
the same total length {Lj), applied electric field {V It ^), and the piezoelectric material constant 

(c/j,) as this nominal actuator, we are interested in developing an unimorph actuator that will 
maximize the value of the non-dimensional dc response, W , while maintaining a minimum value 
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of the natural frequency at 137 Hz. This can be achieved by optimizing the piezoelectric 
thickness (r^), the shim thickness (r^), and the length of the piezoelectric patch (L^). We are 

also interested in estimating the gain in W that can be achieved by developing an optimized 
bimorph actuator that has the same total length, piezoelectric constant and the applied electric 
field as the nominal design, at the same frequency. Therefore an actuator optimization scheme 
based on Matlab® was implemented to maximize the non-dimensional dc response, W. The 
‘finincon’ fimction in Matlab® minimizes the cost function subject to a set of constraints on the 
actuator design parameters, thus giving the maximum attainable value for the cost function. The 
negative of non-dimensional tip deflection, -W was chosen to be the cost function for the 
optimization. The cost function is expressed as. 


3E{l+7)t[L-2)L^ 




{\+£>ET+6EV^M^E+E^t*y 


Similarly, the cost function for the bimorph can be expressed as 


(33) 


6E{\+T)T(L-2)L^ 




(34) 


(l + 6Et +12Et^+8i^E)' 

To gain more insight into the choice of optimization function, let us examine the terms on the 

LHS of Eq. (34) used to non-dimensionalize the tip deflection, . The piezoelectric constant 


djj is a constant for a given piezoelectric material. Similarly the total length of the piezoelectric 
actuators L is also a constant. That leaves the electric field, F/t^, where F is the applied 


external voltage and is the thickness of the piezoelectric patch. It can be seen easily that 
maximiun deflection can be achieved when the value of the electric field equals the value of the 
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coerecive field of the piezoelectric. Thus the choice of the optimization function fV or w , is 

up 

not crucial. 

The design variables involved in the optimization are the thickness of the piezoelectric 

and shim ( ?, ) as well as the length of the piezoelectric (L^). An optimum design for an actuator 

that has a natural frequency value that is at least 137 Hz, is sought. This constraint can be 
represented in the mathematical form as 

f„>UlHz (35) 

The maximum patch length of the piezoelectric should not exceed the total length of the actuator 

and the design parameters should always be non-negative, giving rise to the remaining inequality 

constraints, 

L^<L^ (36) 

X^>0, /^>0, t;>0 (37) 

The optimization scheme seeks the values of , and that maximizes the value of W . The 

result of the optimization study along with the result of the analytical static analysis of the 
nominal design is summarized in Table 2. 
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Table 2: Results of the optimization study. 


Parameter 

Optimized Unimorph 

Optimized Bimorph 

Nominal design 

Lr 

70 mm 

70 mm 

70 mm 

ts 

309 |xm 

98.5 pm 

508 pm 


404 pm 

305 pm 

254 pm 

h 

44.5 mm 

46 mm 

35 mm 

/„ 

137 Hz 

137 Hz 

137 Hz 

W 

63.1 

127.9 

44.2 


It can be seen from the table that optimizing the design variables can increase fV 1.5 times 
the nominal design. The value of W obtained from the optimized bimorph case is almost the 
double of that for the optimized unimorph. In order to validate the choice of the optimization 
function chosen, a separate optimization for a unimorph actuator was done where the cost 
function was instead of W , with the additional constraint that V lt^=V^, where is the 

coercive field and has a value of 30 V/mil. The resulting tip deflection had a value of 
0.00142 m which is the same value for obtained from optimizing W . 

Thus, by constructing an actuator based on the optimized dimensions, maximum tip 
deflection can be attained for a prescribed natural frequency. This result is of particular 
importance in flow control applications, which makes use of large volume displacements of the 
actuator, while placing constraints of the bandwidth of generation. 
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3.8 Conclusions 

An analytical static analysis, and two Finite Element Models, one using perfectly bonded 
shim and piezoelectric elements and the other using a shear element for the bond between the 
shim and the piezoelectric were developed and compared to tip displacement frequency response 
measurements of piezoelectric unimorph flaps. A parametric study was conducted in which the 
size of the piezoelectric patch was systematically varied from 10% to 100% of the entire beam 
length in approximately ten equal steps (while all other variables were fixed). The purpose of 
the study is to develop and validate a design tool, which will ultimately be used to design such 
actuators. 

The dc response as well as the natural frequency of the actuators were extracted from the 
frequency response fimction plots and compared to that obtained from the models. All models 
predicted the dc response accurately. The analytical natural frequency served as an upper limit 
for the theoretical predictions, the best estimate of which was obtained by the shear lag model. 

The benefit of the analytical model is that it describes the fimctional dependence of the non- 
dimensional dc response on design parameters like LJLj, E^IE^, tJLj. and the 

dependence of non-dimensional natural frequency on design parameters like 

Ep! E^, and for both unimorph and bimorph actuators. The analytical model is much 

simpler than both the finite element models and it predicts the response to a reasonable accuracy 
of better than 15%. The analytical model was also used to formulate an optimization procedure 
that maximized the non-dimensional dc response, for a given system bandwidth. 

An in-depth understanding of the structural dynamic characteristics of piezoelectric actuators 
will enable us to reliably utilize them in the design of actuators for flow-control applications. In 
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addition, the important question of coupling between the actuator and fluid motions can be 

rigorously studied. 
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Appendix A - FEM Formulation for the Shear Lag Model 


A schematic of the model is shown in Figure A-1. The flap is discretized into elements in 
the gap region, elements in the composite region (i.e., one each for the shim, bond, and piezo 
portions), and elements in the shim region. 



Figure A-1: Node-numbering scheme for the finite-element, shear lag model. 

The piezo and shim elements are similar to the “Hermitian” beam elements used in the 
perfect bond FEM [Cattafesta et al. 2000], having 2 nodes and 3 degrees of freedom {u,w,6) at 
each node. Here, u is the horizontal or x displacement of the nodes, w is the vertical or z 
displacement, and 6 is the angular displacement. Each element in the bond layer has four nodes 
with two degrees of freedom (w,w) at each node. In the perfect bond model, the nodes are 

defined at the midpoint of the element faces However as shown in Figure A-1, for 

the shear lag model, the nodes need to be offset downwards (upwards) for the shim (piezo) 
elements to allow these nodes to coincide witii the four nodes of the bond element . A 

schematic of the element nodes, degrees of freedom, and transformation is shown in Figure A-2. 
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Figure A-2: Schematic of the discretized flap elements. 
The node offset is achieved by defining a transformation matrix T 


T = 


1 

0 

+1 

2 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

2 
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0 

0 

0 

1 

0 
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0 

0 

0 
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1 


(A.1) 


where the “+” and signs are used to shift the shim and piezo elements, respectively, and t is 
the thickness of the element. 

The elemental mass matrix is derived from the equation for kinetic energy of an element. 
For shim and piezo elements, the transformed elemental mass matrix is given by 


m^=T^fhT 


(A.2) 
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’140 0 0 70 0 0 

0 156 22/ 0 54 -13/ 

_ m/ 0 22/ 4/" 0 13/ -3/' 

m = 

420 70 0 0 140 0 0 

0 54 13/ 0 156 -22/ 

0 -13/ -3/" 0 -111 Af 

where the superscript “T” denotes “transpose,” m is the mass per unit length, and / is the 
length of the element. The mass of an individual bond element is lumped into equal masses at 
the four nodes of the element, giving rise to an 8x8 diagonal matrix with m//4 along the 
diagonal. 

The element stif&ess matrix is determined from the strain-energy relation. For the piezo and 
shim elements, it has the form 



k^ = T^kT 



0 0 0 0 
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12 

6/ 

0 
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6/ 
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6/ 

IP 

0 
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4/' 


(A.4) 


(A.5) 


where E = or E ^ ) 


is the elastic modulus, A is the area, and I is the moment of inertia 


about the local neutral axis. The elemental stif&iess matrix for the bond layer is given by 
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where G and are the bond shear and elastic modulus, respectively, / is the lengUi, b is the 
width, and is the thickness of the bond layer. 

There are no external loads of interest acting on the shim and bond regions. However, for the 
piezo element, the load vector is given by 


F^=rF 


F = 


0 
0 


0 

0 


(A.7) 


(A.8) 


As usual, the elemental mass and stiffiiess matrices are then combined to obtain the global 
mass and stiffiiess matrices M and AT, and the boundary conditions are imposed. To include 
the effects of damping, we use Rayleigh proportional damping (see Thomson, p. 191 for details). 

The frequency response of the actuator may then be calculated by assuming a solution of the 
form Y = subject to harmonic forcing F = F^e““' , where [l^] is the nodal displacement 
vector. The governing matrix equation becomes 
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MY +CY+KY =F 


(A.9) 


The solution to the undamped, free vibration problem is an eigenvalue problem 

[a:-g?^m]To =0 (A.10) 

in which the natural frequencies of the beam are defined by the square root of the eigenvalues of 

[K-a>^M] (A.11) 

The solution to the damped, forced vibration problem is given by the matrix solution 

Yq=]JC+ i(oC - J * Fq (A. 12) 

where co is varied from zero to a specified frequency in discrete steps. 
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